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I.  INTRODUCTION 


The  Fourier  transform  of  a function  I(t)  is  given  by 


S(f) 


00 

-/ 


I(t)  exp  j(-2jtft)  dt 


(1) 


S(f)  will  have  S-function  values  if  I(t)  is  sinusoidal.  Except  for  the 
finite  amplitude,  an  equivalent  expression  can  be  obtained  by  integra- 
tion over  only  one  period  of  I(t),  provided  only  the  values  of  f for 
which  f = m/T  (m  = integer)  are  used: 


I 


S(m  Af)  = / I(n  At)  exp  j(-2nm  Af  t)  dt 


(2) 


where  Af  = 1/T.  Equation  (2)  is  the  basis  of  the  discrete  Fourier 
transform  (DFT) . The  HP  5451B  Fourier  analyzer  evaluates  Equation  (2) 
with  the  integration  replaced  by  a summation: 


N-l 

S ' (m  Af)  = At  y I(n  At)  exp  j(-2jtm  Af  n At) 
n=0 


(3) 


In  the  last  equation,  the  period  is  divided  into  equal  segments 
At  = T/N.  Therefore,  At  must  be  sufficiently  small  so  that  it  will 
resolve  pertinent  structure  (i.e.,  N must  be  large  enough).  Only  for 
this  condition  will  a meaningful  frequency  spectrum  be  displayed. 
Also,  T must  be  sufficiently  large  so  that  the  frequency  resolution 
Af  = 1/T  is  adequate. 


The  Fourier  analyzer  output  F is  actually  1/T  times  S ' (m  Af ) . 

O’ 


This  insures  that  if  an  input  amplitude  at  frequency  n is  In,  then  F 


has  an  amplitude  output  of  Iq/2. 

When  analyzing  a segment  of  a periodic  waveform,  care  must  be 
taken  to  insure  that  the  "window"  is  observed  over  a complete  period 
or  integral  number  of  periods.  If  this  is  not  done,  error  will  result. 
In  the  case  of  far-field  speckle  data,  there  is  no  exact  periodicity. 
The  theory  developed  in  Report  No.  RR-TR-76-2  takes  the  error  induced 
by  integration  over  a finite  window  into  consideration.  The  Fourier 
analyzer  can  be  used  directly  without  concern  about  the  exact  value 
for  the  window.  More  discussion  on  this  matter  is  taken  up  later. 


In  Report  No.  RR-TR-76-2,  the  following  transform  was  defined: 

r1 

F*  (u)  - / I'(x)  [exp  j(-ux)]  dx  , (4) 

-Xl’Xi  J 

1 


3 


where  x = sin  9 (0  = angle  of  detection  with  respect  to  object  surface 
normal).  For  I'(x)  = I'(-x)  the  relation  between  the  real  part  of 
F'  and  the  real  part  of  the  transform  with  integration  limits  from 

-*1’X1 
0 to  x^  is 


2 


Re 


The  DFT  obtained  by  the  Fourier  analyzer  is 
N-l 

F(m  Af)  - I(n  At)  exp  j(-2itm  Af  n At) 

n«0 


(5) 


(6) 


The  desired  DFT  is 


N-l 

2Fq  N_1(m  Au)  = 2 Ax  ^ I'(n  Ax)  exp  j(-m  Au  n Ax) 

n=0 


(7) 


Equivalence  between  Equations  (6)  and  (7)  is  obtained  by  letting 


Ax  * At 

Au  = — Af 
W0 


(8) 


Equation  (7)  becomes 

N-l 

2Fq  N l(m  Au)  = 2(0^  ^ I'(nu)QAt)  exp  j(-2jnn  Af  n At)  . (9) 

n=0 


Defining  I'Ci^n  At)  * I(n  At)  yields 


^O.N-l^”  Au>  = 2w0T  F^m  Af> 
hence , 


(10) 


Re 


2“oT 


Re)  F (m  Af ) 1 

» I 


do 
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The  relation  between  u and  the  variable  representing  displacement  A on 
the  illuminated  surface  is  2k  A = u.  It  is  noted  that  u is  dimension- 
less. From  Equation  (8)  it  is  seen  that  the  real  part  of  the  transform 
F1  as  defined  in  Report  No.  RR-TR-76-2  (essentially  an  autocorrelation 
integral  within  a resolution  length  6^)  is  equal  to  2WgT  times  the  real 

part  of  the  DFT  obtained  by  the  Fourier  analyzer  with  time  and  frequency 
related  to  x = sin  0 and  u = 2k  A by  Equation  (8) . 

It  should  be  noted  that  Ax  = w^At  requires  a variable  angular  scan 
rate.  Because  Ax  = cos  0 A0 , then  d0/dt  = Ax  (cos  0 At)  = Uq/cos  0. 
Hence,  the  angular  scan  rate  goes  as  the  reciprocal  of  cos  0. 


II.  RELATION  BETWEEN  p(b)  AND  THE  RECTANGULAR  WINDOW 


In  Report  No.  RR-TR-76-2,  a function 


F'(u) 


A*A[ exp  j(-ux)] 


dx 


(12) 


was  defined,  where  A*A  = I(x)  as  presented  in  Equation  (4)  and  x 
(0  is  arc  angle).  This  could  be  written  (setting  u = 2jtu') 

00 

F ' (2jru')  = J A*A  W(x)  [exp  j(-2xu'x)]  dx  , 

-00 

where 

II  , -xL  < x < xL 
0 , Otherwise  . (14) 


= sin  0 
(13) 


W(x)  is  the  window  function.  It  is 
transform  of  the  rectangular  window 
A*A. 


seen  then  that  F 1 is  the  Fourier 
function  timet  the  field  intensity 


During  the  derivation  of  explicit  forms  for  F',  a function  p(b/2k) 
where  b = 2k(g  - a + 7tu'/k)  was  introduced.  This  function  was  not  purely 
a window  function,  but  incorporates  it  because  it  was  obtained  by 
integrating  over  the  finite  limits  shown  in  Equation  (12).  The  relation 
between  p and  W is  given  as 


W(x)  [exp  -j2k(£  - a)x] 


(1  - 


(15) 
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where  7 denotes  a Fourier  transform.  It  is  seen  that  p also  incorpor- 
ates a phase  obtained  from  factoring  out  part  of  the  product 

2 2 

A(!,n)A*(a,B) : furthermore,  the  factor  (1  - x ) = cos  9 is 
incorporated. 

For  the  case  of  x^  very  small  and  (■  = a,  p(b/2k)  should  become 

essentially  the  transform  of  the  rectangular  window  which  is  a sine 
function.  It  is  recalled  that 


P(b/2k)  = | 


(16) 


The  significance  of  | = a will  be  taken  up  shortly.  When  x^  -»  0 
and  b is  small.  Equation  (16)  becomes 


This  is  the  same  as  the  limit  of  2x^  sin  bx^/ (bx^)  (2x^  times  a sine 
function)  for  sufficiently  small  bx^.  Now  it  is  assumed  x^  -»  0 and 
b 2:  10.  Equation  (16)  becomes 

2 f ^X1 

p(b/2k)  -»  ( — — cos  bx  + sin  bx^ 


Because  x^/b  is  so  small,  this  further  goes  to 


2 sin  bx. 
p(b/2k)  -* r 


Again  this  is  2x^  times  a sine  function. 

Bie  significance  of  | = a.  can  be  visualized  when  the  distant 
source  is  a point.  For  this  case 


A(pl>  = *0  6(|  - £q)  8(t)  - tiq) 

Now  A(0)  is  given  by 

A(fl)  “ f A(p1)  (exp  jkr)  cos  9 ds 


(17) 
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1 


or 


A(0)  - exp  jkR  / aQ  6(6  - 6Q)  5(rj  - t^) 


X exp  (-jk£  sin  0)  d|  dr)  , 


Therefore,  A*A  is  given  by 


2 

cos  0 


A*A(0)  = 2 j f a0  5^a 

(KR)  trt  a,p 

X exp  [ - jk(6  - a)  sin  0]  d£  d-q  da  d0  . 


!0>  &<e  - n0) 


The  6-functions  force  the  following  result: 


A*A  (0  ) 


Sq  cos  0 


Equation  (19)  was  obtained  because  the  integral  gives  no  contribu- 
tion for  6*0!.  Hence  the  exponential  term  in  Equation  (18)  vanishes. 
Because  p incorporated  this  factor,  it  can  be  neglected  in  Equation  (15), 
For  sin  0 = x small,  p(b/2k)  -*  7 (W(x)}.  Thus,  for  a point  source  and 
observation  window  perpendicular  to  the  distance  to  the  source  (as  well 
as  small)  the  only  significant  values  for  p(b/2k)  become  the  Fourier 
transform  of  the  rectangular  window  function.  As  a matter  of  fact,  F' 
itself  becomes  the  same,  except  for  a factor: 


X1  a0 

f j (exp  - j 2jtu ' x)  dx 

J (\R} 


2 » 

M 


W(x)  (exp  -j2jtu'x)  dx 


F'  (2jtu') 


aQ  sin  2jtu'x1 

7 i point  source  and  x.  small 

\R)2  2irU  X1  1 


It  may  be  of  some  question  as  to  what  effect  the  small  side  oscilla- 
tions of  p(b/2k)  may  have  on  F'  in  general.  The  effect  should  be  to 
decrease  the  value  of  the  correlation  peak  by  a small  factor  and 
decrease  the  shoulders  of  the  correlation  curve  by  a smaller  factor. 

This  amounts  to  only  a small  loss  of  resolution  in  the  correlation 
value  obtained.  Oscillations  will  not  appear  explicitly  because  they 
are  integrated  over. 

Ml.  CONCLUSIONS 


Applicability  of  using  the  Fourier  analyzer  in  far- field 
speckle  measurements  to  obtain  discrete  Fourier  transforms  which  are 
essentially  autocorrelation  functions  with  a resolution  [ dependent 

upon  the  function  p(b/2k)]  has  been  confirmed.  The  manner  in  which  the 
function  p is  related  to  the  transform  window  has  been  clarified.  The 
function  p changes  with  increasing  window  size  to  give  better  resolution 
in  the  desired  autocorrelation  lengths.  The  position  variable  x = sin  Q 
which  is  related  to  the  time  variable  that  is  read  into  the  Fourier 
analyzer  must  change  at  constant  rate.  This  implies  a variable  angular 
scan  rate. 
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